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Abstract 

This paper considers a Josephson Junction array with the geometry of a 
ladder and anisotropy in the Josephson couplings. The ground state problem 
for the ladder corresponds to the one for the one-dimensional chiral XY model 
in a two-fold anisotropy field, which allows of a rigorous characterization of 
the ground state phase diagram and the relevant elementary excitations for 
the system. The approach to equilibrium, which we study using Langevin 
dynamics, shows slow relaxation, typical of systems whose energy landscape in 
the configuration space consists of a wealthy of metastable states, dynamically 
disconnected. 
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I. INTRODUCTION 



Arrays of Josephson junctions (JJ) in the presence of an external magnetic field, are 
probably one of the best examples of physical systems in which the ideas of competition 
between interactions (frustration), disorder effects and complex (glassy) dynamics can be 
theoretically and experimentally checked in a controlled way0. Most of the effort devoted 
in the last decade to these systems has been addressed mainly to the study of two and 
three dimensional arrays as models of extreme type II superconductors in connection with 
some problems in granular superconductors (as the high T c materials)!. In this paper we 
will study a simple geometrical configuration of links, a ladder, which nevertheless shows 
an interesting nontrivial behaviour, and present some results concerning the equilibrium 
(ground state) properties as well as the dynamical approach to equilibrium using Langevin 
dynamics. 

For two-dimensional arrays the ground state is in general unknown for arbitrary values of 
the frustration parameter / (which is essentially the external magnetic field in the appropri- 
ate units). Several approaches have been attempted in order to get close to the truly ground 
state. Halseyi proposed a kind of one dimensional solution (Halsey's staircase) which gives 
correct configurations but only for certain values of /. He found that this solution provides 
a highly discontinuous function for the ground state energy versus /, which is in contradic- 
tion with an exact result indebted to Vallat and Becki. Another way is a numerical attack 
of the problem, by using either phase or vortex variables. Examples of such approach are 
the pioneering works of Teitel and Jayaprakashi and more recently, the "editing method" 
of Straley and coworkers!. In all these methods it is assumed that the ground state for 
rational values of /, let say / = p/q, is periodic with periodicity q (or 2q in some cases) and 
the density of vortices is equal to /. However, for the geometry studied here, because of 
the free character of the boundary conditions on the upper and lower branch of the ladder, 
superconducting currents on them are possible, allowing of a ground state vortex density 
different from /. 
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Different studies on superconducting networks with a ladder geometry are found in lit- 
erature. Fink and coworkers and Simonin et al. have studied a ladder of superconducting 
wires in the Ginzburg-Landau approximation! Current structures commensurate with the 
underlying lattice appear as solutions of the Ginzburg-Landau equations. 

Previous to our work, Kardari first and later Granatoi have study a ladder of JJ in the 
presence of a magnetic field and charging effects. Kardar, by doing various approximations 
to the interaction potential connects the JJ ladder and the Frenkel-Kontorova (FK) model 
with the 'dual Coulomb gas'. Granato focused his study in the effect of small deviations from 
a commmensurate field as a function of the charging energy and in the critical behaviour in 
the absence of external field (quantum XY model). 

In the remainder of this Introduction, the specific model for the Josephson Junction lad- 
der will be introduced, along with the notation used and the main approximations leading 
to it. In section II, making a judicious choice of the gauge, it is established the equiva- 
lence, regarding the ground state problem, with the one-dimensional chiral XY model with 
anisotropy. Going further, an important argument due to Griffiths and ChouS and Sasaki 
and GriffithsO, allows the applicability of the vast amount of rigorous results on the ground 
states of models of spatially modulated structures with convex interactions^ to our system. 
This provides in a rigorous way the main properties of the ground state phase diagram, 
some of which were already suggested by Kardar and Granato. We also use the effective 
potentials method0 for the explicit computation, with arbitrary numerical accuracy, of the 
ground state configuration for any values of the model parameters. 

In section III we perform a linear stability analysis of representative ground state struc- 
tures which reveals its persistence as metastable states outside the domain of / values for 
which they are ground states. Many other metastable structures do exist in the system 
and the dynamical approach to equilibrium, which we consider using Langevin dynamics, is 
characterized by a constrained dynamics^ leading to slow relaxation. This is a remarkable 
behaviour in the absence of quenched disorder; i.e., no randomness in the Hamiltonian of 
the system0. 



The system we consider is a ladder of superconducting islands in the presence of a 
magnetic field perpendicular to the plane of the ladder (see Fig. [I]) and we assume that each 
island is proximity coupled to its three nearest neighbours. The interaction Hamiltonian for 
the system is 

H = - Yj Jij c os{di - 6j - A i:j ) = - Jijcos^ij (1) 

<ij> <ij> 

where 9i denotes the phase of the superconducting order parameter at the ith island or 
site; jij, the gauge invariant phase difference, is restricted to the interval (— tt, it]; and is 
proportional to the line integral of the vector potential A between the ith and the jth sites, 

2tt pi -» -» 

Mi = — A-dl (2) 

It is required that J2p Mj = 27r/, where / is the ratio of the flux caused by the external 
field with the superconducting magnetic flux quantum and is a measure of the frustration. 
This relation expresses the discretized Maxwell equation for the vector potential and the sum 
is taken in a clockwise direction over the bonds surrounding the plaquette p of the lattice. 
Because the phases are all defined in the interval (— n, it], we have J2 P lij — 27r(n — /) where 
n is the integer that defines the vorticity in each plaquette. Associated to this value we 
define the vortex density u as the mean value of n in the ladder. 

The hamiltonian (Eq. [I]) is the sum of the Josephson coupling energies between the 
neighbouring islands. Here we are neglecting screening currents by assuming that Aij is 
fully determined by the external magnetic field. This assumption is correct whenever the 
penetration length is much greater than the island width. In this situation, there is no flux 
quantization!!!!, the vortex density being not a flux quanta density, but a fluxoid quanta 
density. Also we consider that no charging effects are present. For the coupling constants, 
Jij, we will assume that J^- = J x for horizontal links and J^- = J y for vertical ones. 
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II. GROUND STATE AND PHASE DIAGRAM 



Among the different choices for the gauge one is particularly convenient: A^j = +ivf 
for the upper links, Aij = —irf for the lower ones and Ay = for vertical links, which 
corresponds to a vector potential parallel to the ladder and taking opposite values on upper 
and lower branches. Thus 

H = - Ei [J x cos(9i - 9 i+1 - Ti f) 

i ^ 

+J x cos(9' i - i+1 + irf) + J y cos(9i - 9\) 

Here 6 l i (6 1 -) denotes the phase on the upper (lower) branch of the ladder at the ith step. 
It is easy to see that the phase configurations which minimize the hamiltonian are such 
that 9i + 9 i = constant, independent of i; then, by fixing this constant to and defining 
the anisotropy parameter, h = J y /2J X , one obtains the equivalence between the following 
hamiltonian, 

H = -2J X ]T [cos^i - 9 i+1 - irf) + hcos{29 i )\ (4) 

i 

and(Eq. ^|), regarding the ground state configurations (and other local minima for the 
energy) . 

The hamiltonian (Eq. [|) describes a one-dimensional chiral XY model in a two-fold 
anisotropy fieldS. It belongs to a general class of one-dimensional models of spatially modu- 
lated structures0, the simplest of them being the FK model, extensively studied by AubryEl. 
The equilibrium properties of these models depend crucially on the convexity of the interac- 
tion potential. In the model defined by Eq. f|, the nearest neighbour interaction potential is 
nonconvex. However, it can be proved00 that only the convex part of the interaction term 
plays a role in determining the ground state configurations, so that the ground state prop- 
erties are those of a convex model; of course, the situation may be different if one considers 
different aspects of the model, other than the ground state. 

The essential physics of the model is the competition between the anisotropy term (com- 
ing from the vertical Josephson couplings) which tends to pin the phases value to or tt, 



and the interaction term (coming from the horizontal Josephson couplings) which tends to 
fix 9i — 9 i+ i = 7t/, that is, it tries to keep the value of the vortex density at the frustration 
value, u = f . The ground state phase configuration at given values of the parameters (h, /) 
is the result of the compromise between both competing tendencies. In order to compute 
them we use the effective potentials methodS which has become the standard method to 
obtain the phase diagram for this type of models. When used in combination with Newton 
and relaxation techniques it gives the ground state configuration with arbitrary accuracy. 
The computed phase diagram is shown in Fig. |2| where, for clarity, only a few transition lines 
are represented. Characterizing a ground state by the value of the vortex density (a;), we 
see in the Fig. ^, as predicted by the rigorous arguments mentioned aboveffl, that u(f), 
for fixed h, is a Devil's staircase: a continuous function but such that for each rational value 
of the vortex density there is an interval of values of the field for which u remains constant. 

This phase diagram is quite different from the one expected for the isotropic (h = 0.5) 
two-dimensional J J arrays. For the 2D system, though there is no rigorous proof for it, it 
is assumed that there are no intervals of stability for rational values of u, that is, oj = f 
everywhere!!!. In the case of the ladder, however, the vortex density is not equal to the 
field. If we move along inside a step, the ground state configuration for the gauge invariant 
phase differences do change with /, while the vortex configuration remains unchanged; 
correspondingly, supercurrents along the links of the ladder keep varying to compensate 
for the increase of the field with no change in the vortex (fluxoid) density. In the case of 
the ground state for zero vortex density, it is tempting to speak of this effect as a sort of 
Meissner effect, but one should notice that no flux is expeled from the ladder, and one cannot 
interpret the value of / for which the ground state changes as an analog of the critical field 
for superconductors: it is the fluxoid density what changes from zero at that value, not the 
flux density. 

Although uj(f) shows such complex aspect the ground state energy is a continuos function 
of the frustration (Fig. |3|). This is also the case for a 2D J J array, as proved analytically in 
Reft 
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The vortex configuration rij corresponding to a ground state of vortex density u(< 1) is 
explicitly given by 

rii = Xuj(iu + a) (5) 

where a is an arbitrary constant and Xw{ x ) = Xu>i x + 1) is the characteristic function of the 
interval [0, to): 

{1 if < x < uj 
(6) 
if uj < x < 1 

Then, rij is a periodic sequence (with minimal period) for rational values of uj and a 
quasiperiodic sequence for irrational values of the vortex density; it is traditional to speak 
of commensurate and incommensurate ground states respectively. 

Commensurate ground states are, for any value of the parameter h, pinned and defectible. 
A configuration is pinned when there exists a finite value Id (depinning current) such that 
if a current I < Id is injected into each island on the upper branch and extracted from each 
island on the lower branch, the vortex configuration remains unchanged. In this case the 
phases change to a new equilibrium configuration and no voltage appears on the links (the 
ladder remaining superconducting). For values of the external current greater than the de- 
pinning current, the phases configuration varies with time and a voltage can be measured. A 
defectible configuration admits discommensurations (defects). An elementary discommensu- 
ration in a commensurate configuration corresponds to a domain wall separating equivalent 
vortex configurations which are shifted relative each other, with minimum increase (or de- 
crease) of the local vortex density (see figure |] for an example). Notice that the vortex 
configuration of an elementary discommensuration, though only locally different from the 
underlying commensurate vortex configuration, cannot be obtained from this through a fi- 
nite number of local changes, but entails the whole rearrangement of a semiinfinite part of 
the system. It is important to keep this point in mind when dynamical approaches to equi- 
librium are studied. The creation energy of an elementary discommensuration goes to zero 
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at the value of the frustration parameter where the ground state vortex density changes, 
and a C-IC transition takes place. 

Incommensurate ground states show two different regimes, separated by an Aubry tran- 
sition (Transition by Breaking of Analiticity)0 at a critical value h c of the parameter h, 
which depends on the irrational vortex density, u. Below this critical value, the ground 
state configuration is unpinned (any external current produces the appearance of voltage on 
the links) and no defects can be sustained. In this regime the sequence of gauge invariant 
phase differences 7^ can be expressed in terms of an analytical hull function. Thus, in the 
case of a vertical link 7$ = 9i — Q i — g(—i2iruj + (3), with (3 an arbitrary constant (Fig. ^|). 
The situation changes when h grows above the critical value h c : the hull function develops 
infinitely many discontinuities and the incommensurate ground state becomes there pinned 
and defectible. Our estimate of h c for a golden incommensurability ratio, to = (3 — y/E)/2, 
is h c = 0.245... . This estimate should certainly be improved, for we have used rather poor 
rational approximants of uj. On the basis of the plausible irrelevance of the deviations from 
quadratic of the interaction potential, one can conjecture that the main gap of the hull 
function for golden irrational vortex density behaves as A ~ {h — h c p, with £ = 0.712, the 
critical exponent obtained by Mackay0 for the Aubry transition in the standard FK model. 
The depinning current, I4, for the golden incommensurate structure has been estimated us- 
ing simulations in the RSJ approximation!^. It has been shown that scales as (h — h c ) v 
with v = 2.75 close to the estimation, v = 3.011, of MacKay for the standard FK model. 



III. METASTABILITY AND RELAXATION PHENOMENA 

One of the characteristics of frustrated models is the existence of a large number of 
metastable states, a feature which influences dramatically the dynamical approach to equi- 
librium. Those states are local minima of the energy (i.e.: stable solutions of the equilibrium 
equation, |^ = 0). For the model we are considering, it seems plausible that the existence 
of truly chaotic metastable states could be justified, as an extension of the results obtained 



8 



in Ref£2l, for the FK model. In order to illustrate the nature of metastability in the JJ lad- 
der, we will consider the linear stability of a very restricted class of configurations, namely 
those vortex configurations which are ground state for some values of the model parameters. 
To analyse the stability of such states we have worked out the spectrum of small linear 
perturbations. The procedure has been the following: 

a) First, by fixing the model parameters at values inside the tongue (see phase dia- 
gram) corresponding to the selected particular value of u>, and using the method of effective 
potentials, the ground state phase (and vortex) configuration is obtained. 

b) Now we vary finely the parameter / (typically Af ~ 10~ 3 ) and using a Newton 
method to solve the system of (nonlinear) equilibrium equations |^ = 0, we determine the 
evolution of the equilibrium vortex configuration under quasistatic changes in /, along with 
the energy variation. 

c) At each value of /, the spectrum of the small perturbations matrix, { qq.qq. }, around 
the corresponding equilibrium phase configuration, is computed. 

In all the cases, a zero eigenvalue is found, which corresponds to the (continuous sym- 
metry) invariance of the Hamiltonian (^) under uniform rotation of all the phases. If the 
rest of the eigenvalues are all positive, the configuration is linearly stable, and the state 
is a local minimum of the energy for the fixed value of u under consideration; when the 
lowest eigenvalue takes on a negative value, the configuration is linearly unstable, usually 
corresponding to a local maximum of the energy. 

In Fig. ^ we show the evolution of the (non-zero) lowest eigenvalue of the stability matrix 
for some simple values of the commensurability ratio u = 1/2, 1/3, 1/4, 1/5,2/5 when / is 
varied between and 0.5 and h = 0.5 (J x = J y ). Not surprisingly, the range of stability is 
much wider than the interval of / values for which a given state is a ground state: a C-IC 
transition does not have associated a lack of stability of the commensurate state (which 
remains a local minimum of the energy), but it corresponds to the vanishing of the creation 
energy of elementary discommensurations, a feature that cannot manifest itself through a 
linear stability analysis. At the edge of the stability intervals, marked by the negative sign of 
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the lowest eigenvalue, the extremal character of the configuration changes from a minimum 
to a maximum of the energy. An interesting feature is that, at the border of the stability 
interval, the gauge invariant phase configuration is always such that there is a link where 
the supercurrent reaches its maximum (critical current) value, which in turn coincides with 
the interaction potential leaving its domain of convexity. At this point any small change in 
the field cannot be sustained by an increase of the currents. The vortex structure becomes 
unstable and the nearer (in configuration space) stable phases possess a different value of uj. 

In Fig. we represent, as a function of the frustration, the energy of some simple com- 
mensurate states. The ground state energy corresponds to the lower envelope of these (and 
infinitely many other) curves. The stability transition points are marked by open circles. 
We can observe that near the borders / = 0or/ = l/2 the energies of these commensurate 
states are well separated while at intermediate values of / they are very close. The energies 
of other (commensurate and incommensurate) stable states, not included in the figure, also 
lie around. Besides all those stable states corresponding to minimum energy configurations 
which exist as ground states for some / values, one can construct other metastable states, 
corresponding to almost arbitrary rearrangement of vortices, with energies also lying about. 
Then, for intermediate values of /, the energy landscape consists of an extremely complex 
set of local minima, with comparable energies, corresponding to phase configurations which 
are, generically, rather separated in configuration space and, from a dynamical perspective, 
almost disconnected, a situation which is sometimes referred to as constrained dynamic^. 
Remind that most of the states are unreachable from a given one. The dynamically reachable 
states are those which come from the annihilation (or creation) of a finite density of vortices 
and not from the rearrangement of part of the lattice. This fact introduces a hierarchy of 
states in the relaxation dynamics which is relevant in the glassy properties of the modefi 
We can conclude that this model presents ingredients of systems with glass behaviour: a 
complex structure of metastables states and constrained dynamics. It is worthwhile to em- 
phasize that here there is no quenched disorder in the hamiltonian. As we will see such 
scenario is strongly confirmed by numerical simulations of the dynamics in the presence of 
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noise. 

The existence of strong dynamical constraints in configuration space is a microscopic 
feature which leads to the macroscopic phenomenon known as slow relaxation. Such be- 
haviour has been observed in many systems like spin glass compounds, polymers, granular 
superconductors, etc. In order to check this phenomenon in the JJ ladder we study the 
Langevin relaxational dynamicsil 

a tt 

kit) = -r— + x t (t) (7) 

where H is defined by Eq. |3| and we use T = 1 and J x — 1; Aj is an additive thermal noise 
in the phases and satisfies (A,(t)) = 0, (Xi(t)Xj(t )) = 2T5ij5(t — t). 

We look for the relaxation of commensurate as well as random vortex configurations for 
different values of the temperature T. For each one we have computed the density of vortices 
as a function of time and followed the time evolution of the corresponding vortex spatial 
configuration. In these simulations we use ladders with 400, 2000 and 5000 plaquettes to 
avoid finite size effects. The results presented here have been obtained for / = 0.25 and 
h = 0.5 (J y = J x ). The ground state for those parameter values has zero vortex density, 
u = 0, and the configuration space shows there an extremely complex structure of metastable 
states with close energy values. 

In the case of an initially ordered configuration, three temperature regimes are found 
(Fig. H). At very low values of temperature the state, of course, remains as a metastable 
configuration. From some higher value of T, which depends on the particular initial state, 
and until T ~ 0.05 one observes the decay from the initial ordered state to a disordered 
metastable vortex state. Such state is basically of the same type obtained after the T = 
relaxation of a random phase configuration. For values of temperature above T = 0.05 the 
states relax slowly to the zero vortex ground state configuration. Such relaxation for the 
values of the parameters chosen is observable in the range of temperature from T ~ 0.05 to 
T ~ 0.15. About this last value of T, thermal activation of vortices and antivortices appears 
superposed to the purely relaxational effects. 
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Many classes of functions have been proposed to fit a slow relaxation curve. One of the 
most general is the KWWtl (Kohlrausch, Williams and Watts) or two-parameter stretched 
exponential law, defined by n(t) ~ exp[—(t/r)P]. Slow relaxation corresponds to values of 
the parameter f3 < 1. For (3 = logarithmic relaxation is seen. A value of (3 from 0.5 to 
0.7 is common in glasses. The fittings of our simulation data give values of (3 between 0.5 
and 0.9, depending on temperature (see Fig. |9|). We have not been able to find a functional 
dependence of the exponent (3 with the temperature. However r ~ exp(a/T) as could be 
expected^. 

Finally, we have also investigated the microscopic characteristics of the vortex dynamics 
in the different regimes of the relaxation. At low temperatures the relaxation of the com- 
mensurates phases is dominated by the nucleation of new structures compatible with the 
initial one. The long time state is a metastable one formed by different commensurate struc- 
tures separated by domain walls, see Fig. |H]. The metastable states reached from random 
phase configurations are essentially the same. At higher temperatures these intermediate 
metastable states decay to the ground state configuration. As a result of the characteristics 
of configuration space in this range of temperature we find slow relaxation. Vortices are 
slowly (thermally) expelled from the array, see Fig. [XT]. 

At very high temperatures the dynamics is dominated by the thermal generation of vor- 
tices and antivortices. Phases are randomly distributed between and 2tt and n approaches 
to /. 

IV. SUMMARY AND CONCLUDING REMARKS 

In this paper we have analized a Josephson Junction ladder in the presence of a perpen- 
dicular magnetic field. The ground state problem of this system is equivalent to the one of 
a FK model with convex interparticle interaction, which allows to apply the Aubry theory 
for this class of models. We have calculated exactly the ground state phase diagram which 
shows tongues of stability for rational values of the vortex density and a devil's staircase 
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structure. Incommensurate structures exist and can be described by a hull function. Below 
a critical value h c the hull function is continuous and the structure is sliding (zero depinning 
current). Above such value the hull function develops discontinuities and a non zero pinning 
force appears. Modern microlithographic techniques along with methods to detect vortices^ 
may check the above results in large JJ ladders. 

Relaxation of arbitrary vortex configurations fits to a "slow dynamics" function although 
no structural disorder is present in the model. Disorder is introduced via initial random 
configuration and thermal fluctuations. The existence of a complex structure of multiple 
metastable states and a strongly constrained dynamics in the configuration space are the 
essential ingredients for this "glassy" dynamics. 

The results summarized above have several consequences. Due to commensurability 
effects, finite size in the ladder direction can produce changes in the dynamical response 
under dc and/or ac currents^. Mismatch in the boundary condition generates defective 
vortex configurations with a peculiar dynamics^. This is an important effect in order to 
interpret correctly dynamical results in the ladder. Though the equilibrium properties of 
the JJ ladder are those of a FK model with convex interaction, this equivalence does not 
hold when dealing with the dynamics of the model. For instance, the dc driven dynamics 
in the ladder does not, in general, shows a unique V(I) as it should occur for a convex 
modelil. Finally, if screening currents are considered, the anisotropic JJ ladder can be a 
good model for long Josephson juntions and stacked JJ0. In this way, an extensive study 
of the Josephson Junction dynamics under dc + ac driving currents is in progress. 
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FIGURES 

FIG. 1. Schematic picture of the J J ladder. The gauge choice is shown in the right-most 
plaquette. Eq. ||| gives the interaction hamiltonian of this system. 

FIG. 2. The phase diagram of the J J ladder obtained using the method of effective potentials. 
Each phase is defined by the value of u and, for clarity, only a few of the transition lines are rep- 
resented. The sketchs of the ladders show the vortex configurations of the simplest commensurate 
states. 

FIG. 3. u)(f) for the ground states configurations when h = 0.2. This function is a devil's 
staircase: a continuous function with an step for each commensurate value of uj. The inset shows 
the continuity of the ground state energy as a function of the frustration. 

FIG. 4. (a) An elementary discommensuration (DC) in a uj = state. We show the DC from 
both phase and vortex points of view, (bl) The two possible vortex sequences for the uj = 1/2 
ground state. (b2) Elementary DC in the uj = 1/2 commensurate configuration. 

FIG. 5. 2-7T modulo representation of the hull function for the gauge invariant difference of 
phase in the vertical links for the co = 13/34 ground state. This state reflects the behaviour of true 
incommensurate (irrational uj) phases. Above a critical value of h the function develops infinite 
discontinuities being analitic for values of h lower than the critical one. In this case h c — 0.24. 

FIG. 6. Lowest of the non zero eigenvalues of the matrix of stability of a state when / is varied. 
The picture shows the cases of some simple commensurate states when h = 0.5. Open circles show 
transition from a minimun to a maximun of the energy when / is decreased. 

FIG. 7. Energy diagram of some commensurate states as a function of /. The lowest envelope 
is the ground state energy. Open circles mark the limit of stability of each state, see Fig. |6|. The 
dot line is for / = 0.25, the value of the field we choose in the relaxation calculations we present 
in the Section III. h = 0.5. 
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FIG. 8. Temporal evolution of the mean number of vortices in the ladder, n, in the relaxation 
curves at different temperatures. The initial state is the u = 1/2 ordered state. This state is stable 
at this value of / and h (/ = 0.25, h = 0.5). At low values of T no relaxation is seen. For values 
of T below T ~ 0.05 the state decays to a new metastable one. At temperatures above this value 
it decays slowly to the uj = ground state. At highest temperatures the thermal generation of 
vortices is dominant, being < n)(t) >t= /• In the cases of T = 0.1 and T = 0.15 we also show 
the relaxation curves of random initial configuration. Those curves are nearly close to the uj = 1/2 
relaxation curves. 

FIG. 9. Slow relaxation curves with random intial conditions have been adjusted using the 
two-parameter stretched exponential law in the range of temperatures between 0.05 and 0.15. 
Different exponents between 0.5 and 0.9 are found, (h = 0.5, / = 0.25) 

FIG. 10. Temporal evolution of a vortex configuration in the ladder, rij(t), showing nucleation 
processes. Such processes dominate the 'low temperature' relaxations of ordered structures. Here 
the n{t = 0) = 1/2 ordered state decays to a disordered metastable state, see Fig. ||. In the picture 
a black mark represents a vortex, (h = 0.5, / = 0.25, T = 0.01) 

FIG. 11. Temporal evolution of a vortex configuration in the ladder, rij(t), showing slow re- 
laxation. Such processes domine the 'mean temperature' relaxations. Here the initial state is a 
n(t = 0) = 1/2 ordered state which decays quickly to a disordered metastable state, see Fig. ||, in 
the manner shown in Fig. [l(]. Then, such states decays slowly towards the corresponding ground 
state, in this case (h = 0.5, / = 0.25, T = 0.1) and uj gs = 0, no one vortex in the ladder. 
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